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Using conformal field theory methods Eggert and AmeclJ3 have shown that the semi-infinite 
S = 1/2 Heisenberg antiferromagnetic chain exhibits a remarkable alternation in its local response 
to a uniform field at low temperatures. Such alternation is not an essentially quantum effect: similar, 
and sometimes stronger, susceptibility alternation is a feature of classical Heisenberg-Ising chains at 
T = 0. In S = 1/2 chains, susceptibility alternation is not unique to the Heisenberg model, but can 
be seen in expansions about the Ising model and also in the XY model. 



I. INTRODUCTION 

In a system of coupled spins which lacks translational 
invariance, the response to a uniform magnetic field also 
lacks translational invariance. Labeling the spins by a 
site index n we can define zero-field local susceptibilities 



X n 



a/3 _ 



8(S°) 



(where a and (3 denote Cartesian axes) which will be the 
principal subject of this paper. 

Eggert and AffleckEl recently considered the local sus- 
ceptibilities in the semi-infinite S — 1/2 Heisenberg an- 
tiferromagnet 

n = J2 Js n -s n+1 + h-s„ (2) 

using powerful methods of boundary conformal field the- 
ory. In this case spin-rotational symmetry implies that 
X'/f = 3u,f3Xn- Their results have the following form, 
which we expect to be asymptotically exact in the large- 
n and low-temperature limit (up to possible logarithmic 
corrections) . The local susceptibility can be decomposed 
as Xn = X™ u + ( — l)" _1 Xn lt where x um and x alt are sepa- 
rately smooth, positive functions of n. The "uniform" 
part approaches the susceptibility of an infinite chain 
with increasing n. The "alternating" part is given by 

xf oc n ((3 J 3 sinh(4n//3J))~ 1/2 (3) 

with [3 — ksT and a proportionality constant of order 
unity. This result has the remarkable feature that at 
T = the alternating part increases with n, and indeed 
it diverges like n 1 / 2 . 

The local susceptibility in spin chains is a question 
of more than purely theoretical interest. It is the local 
susceptibility which is relevant to local magnetic probes 
such as NMR, and in any real sample of a quasi-one- 
dimensional magnet there will be some defects that break 
the chains. Indeed, the effects of open chain ends may 
have already been seen in NMR studies of S^CuOaJ 
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The questions we would like to address are complemen- 
tary to those dealt with by Eggert and Affleck. Is there 
a simple way to understand the remarkable local suscep- 
tibility alternation which they found? Is that alternation 
unique to Heisenberg models, or is it a ubiquitous fea- 
ture of semi-infinite antiferromagnetic spin-chains? We 
will explore these questions at T = 0, where the alterna- 
tion phenomenon is strongest. 



II. CLASSICAL SPIN-CHAINS 

Eggert and Affleck suggested that the local suscepti- 
bility alternation was a purely quantum mechanical ef- 
fect. Here we will argue that, at least in the context of 
spin operators, it is actually a classical effect. Indeed, 
we will show that classical version of the model (g) has 
even stronger susceptibility alternation than the 5=1/2 
quantum result (Q), and thus quantum fluctuations, like 
thermal fluctuations, act to suppress the susceptibility 
alternation. 

Let us consider an TV-site classical Heisenberg-Ising 
chain 

AT-l N 

H = J Z S*S* +1 + Jj_(S*S* +1 + S^S^ +1 ) + h • ^ S„ 

n— 1 n— 1 

(4) 

where the spins are unit vectors. We shall ultimately be 
interested in the N — > oo limit. It will prove convenient 
to represent the spins in terms of their azimuthal angle 
<t> and projection on the xy plane s = sin0, and to define 

A = J±/J z . (5) 

If we are going to consider the T = local suscep- 
tibilities it should be evident that for Ising anisotropy 
we should be examining the xx (or yy) susceptibilities 
while for XY anisotropy the zz susceptibilities are of 
interest. In the limit that A = 1 the choice of suscep- 
tibility is irrelevant (so long as it is diagonal) but for 
A < 1 the zz susceptibility is very different in character. 
For even N there are two degenerate ground states in 
zero field (the "Ising states" in which the spins are per- 
fectly aligned with the z axis and antiferromagnetically 
ordered) which remain the ground states until a finite 
value of h z is achieved, so the zero-field local susceptibil- 
ities are all strictly vanishing whether one averages over 
both ground state configurations or not in carrying out 
the T = average. For odd N the two Ising state ener- 
gies vary linearly and cross when h z = 0, so the natural 
characterization of the small-field behavior of Xn would 
be (—l) n ~ 1 25(h) and the zero-field limit would be ill- 
defined. Clearly the N — > oo limit is also ill-defined at 
T = 0. For A > 1 the xx and yy susceptibilities are anal- 
ogously pathological. In contrast, the zero-temperature 
local transverse susceptibilities \% x at A < 1 are well de- 
fined for all finite-length chains and have smooth N — > oo 
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limits; and likewise for Xn Z at A > 1. Let us now calcu- 
late them for general A and N, beginning with the case 
of Ising anisotropy. 

For h = hx the ground state has alternating S z so 
that calculating the ground state configuration involves 
minimizing 

N-l 



w =E(- J ^\/( 1 - s «)( 1 - s «+i) 

n=l 

N 

+ J±S n S n+1 cos(</>„ - <j>n+l) ) -/l^S„COS0„ . (6) 



71=1 



Setting dH/d(j) n — dH/ds„ — leads to the conditions 

- hs n wa.<f> n /J± 

= Sn-iSn sin(0„_i - 4> n ) + s n s n+ i sin(0„ - </>„+i) (7) 
and 



v 1 _ S n-1 + V 1 ~ S n+1 

hcos(j) n /J z = S n 

v 1 - 4 

+ A(s„_i COs(0„_i - </>„) + S n+ i cos(0„ - <j> n +l)) (8) 

for n = 2 . . . iV — 1. For n = 1 (iV) the terms containing 
so and </>o (s^v+i an d <^JV+i) are dropped but the equa- 
tions are otherwise identical. The equations (0) have the 
solution tp n — for all n. This makes perfect sense: why 
should the spins deviate from the xz plane? (Note that 
4> n = for all n will only correspond to the minimum en- 
ergy configuration when h > 0; for h < 0, <fi n = ir gives 
the minimum. There are also solutions to (0) in which 
some of the cf> n are and the others are ir, but these are 
not energy minima.) Then in (|^) all of the cosines can 
be replaced by 1, and s = S x . To obtain the zero-field 
local susceptibilities we linearize (||) with respect to the 
s n and obtain (with x = X zz ) 

2Xn + A(Xn-l +Xn+l) = 1/Jz (9) 

for n = 2 . . . N — 1, and the boundary conditions 

Xi+Ax2=l/J 2 Xn + Axat-i = l/Jz ■ (10) 

Before we proceed to discuss the solution of the lin- 
ear system which governs the \n let us consider what 
happens for XY anisotropy and fields h = hz. In 
this case the ground state will have non-alternating S z 
so instead of (^|) the quantity to minimize will have a 
plus rather than minus sign associated with the term 

! J(1 — sLj). In addition, the external field 



term hs n cos <j) n is replaced by hy/l — s^. Then the 
energy-minimizing solution to the analog of Eq. (0) is 
(j) n — 4> n +i = Ti"- There is clearly degeneracy in the ground 
state: starting from any spin configuration, all spins can 
be rotated by an equal amount about the z axis without 
changing the energy. Then, defining c = cos 9 — S z , the 
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analog of Eq. (@) becomes (extremalizing the energy with 
respect to the c„ rather than the s n ) 

\h _ c n-l + \A ~~ C n+1 
h/J z = Cn-l + Cn+1 - Ac„ 

(11) 

which upon linearization with respect to the c n leads to 
the relation 

2AX„+Xn-1 +Xn+1 = l/Jz (12) 

where x — X zz here. This has the same form as Eq. (0), 
if in the former one makes the substitutions A — ► 1/A 
and J — > AJ. In other words, the equations are related 
by the interchange J z <-> Jj_. So, solving Eq. (JoJ) for 
the transverse local susceptibilities in Ising-anisotropic 
systems is enough: there is no new calculation to be done 
for XF-anisotropic systems. 

Now let us get on with the business of solving for the 
local susceptibilities. Two special cases are of particular 
interest. 

For the Ising model, A = 0, one has Xn — V(2«7*) for 
all interior sites on the chain, while xi = Xn = 1/ 'Jz- 
Note there is no alternation in this case, and for the in- 
terior sites the response is the same as if there were no 
boundaries at all. 

For the Heisenberg model, A = 1 and J z = J± = J, 
the behavior of the system depends markedly on the par- 
ity of N. If N is odd, then the ground-state energy 
varies linearly in \h\ because the ground-state configura- 
tion consists of aligning the odd-numbered sites with the 
field and anti-aligning the even-numbered sites. There- 
fore the local susceptibility has <5-function character at 
every site, with amplitude proportional to (— This 
corresponds to the greatest conceivable local susceptibil- 
ity alternation. For even N, the spins instead cant in 
response to a field, and the local susceptibilities are well 
defined. However, the linear system (^) and ( |l0| ) does not 
have a unique solution: the linear operator has a zero- 
mode of the form (1, —1, 1, — 1, . . . , 1, —1). This apparent 
degeneracy of spin configurations in the small-/i limit is 
resolved by the fourth-order terms in the Si in the Hamil- 
tonian. One finds then that the Xn satisfy Xjv+i_ n = Xn 
and are conveniently described by the recurrence relation 
JXn+i — JXn + {— — n) with the starting value 
Jxi = N/4. For example, for N = 8 the sequence of 
JXn is (2,-1,1,0,0,1,-1,2). There is certainly susceptibil- 
ity alternation here, but unlike the semi-infinite quantum 
chain the alternation is largest near the boundaries rather 
than in the interior. Despite the very different character- 
istics of the local susceptibility for odd and even N, it 
seems that the N — > oo limit for Xn at fixed n is well 
described by (— ) ,l_1 cx): a striking result! 

Let us now turn to the general problem of calculat- 
ing Xn X for < A < 1 and general N, with particular 
attention to the N — > oo limit. We have employed a 
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generating-function technique, which involves extending 
the n values under consideration from zero to infinity, 
defining X(z) — J2^>o X« z ~™- This sort of generating 
function is conventionally known as a "z-transform" .1 In 
order to satisfy (|9|) when n = 1 as well as the lower 
boundary condition in (10) the relation xo — — Xi/A 
must hold; the value of xo (and thence all the \n) will be 
uniquely determined later, in order to satisfy the upper 
boundary condition. Then we can apply the standard 
result for z-transform solution of linear recurrences: 



X(z) 



R{z) + I{z) 
Q(z) 



(13) 



where Q(z) encodes the homogeneous terms in the recur- 
rence (and equals z 2 + 2A~ 1 z+l in the present case), R(z) 
encodes the inhomogeneous terms (and equals [J 2 A(1 — 
z)] -1 here) and I(z) encodes the initial conditions (and 
equals xo^ 2 + (xi + 2A~ 1 Xo)2 = XoQ - Xo + Xi z here). 
In order to carry out the back-transform for X it is 
useful to introduce a partial-fraction decomposition of 
1/Q, which makes it useful to denote the roots of Q by 
z± = A _1 (— 1 ± yl — A 2 ). Then straightforward alge- 
bra, combined with the knowledge that the z-transform 
of (a n — b n ) /(a — b) is z/(z — a)(z — b) and also the relation 
between \o and \i given above, leads to 



1 



AJ z (z + - Z-) 



1 - Z T j 



1 



(14) 



Inserting this expression into the upper boundary condi- 
tion yields the explicit expression for xa 



Xo 



-4 V ~ 2 (A + (1 - A)z+ - z 2 ) + zi v ~ 2 (A + (1 - A)z_ - z 2 ) + (1 + A)(z 4 



2J Z (1 + A) 2 (l-A)(z™ 



JV-l 



(15) 



In taking the N — > oo limit we apply the fact that 
\z + \ < 1 and |z_| > 1, giving 



lim Yn = — ttt - 

N^oo M> A 2 (l + A) 



1 1 

A ~ 2 



v 7 ! - A 2 (\ 1 
1 + A { A + 2 



1 



(16) 



In conjunction with Eq. (|1J), this gives an explicit solu- 
tion for Xn for the semi-infinite chain. 

The solution at hand is explicit but rather unwieldy. 
It turns out that for A not too close to 1 and n not too 
large it is more convenient to apply numerical methods 
directly to the linear system for the Xn- Numerical so- 
lution is especially easy in this case because the problem 
is a tridiagonal linear system. In Fig. [I] we present re- 
sults for xi through X6', although the calculations were 
done for TV = 100, on the scale of this plot any finite-iV 
effects are invisible. Also shown, as the solid line, is the 
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susceptibility of a uniform chain. Some key features of 
the results for A < 1 should be evident: the amount of 
alternation decreases as one moves towards the chain in- 
terior, with the rate of decrease depending on A (smaller 
A giving a stronger decrease), and as A — * 1 all of the 
Xn diverge. What cannot be easily seen in the numerical 
solution are the details of the behavior near A = 1, since 
the finite N effects are strongest there. For this reason it 
is informative to construct an expansion for Xn in pow- 
ers of 8 = 1 — A, using the explicit solution. The leading 
terms are found to be 

Xn « (_)n-i 2 -3/2j-i/2 + (i + 2 (n - l)(-)")/4 (17) 

with the next term OiS 1 / 2 ). Note that up to the order 
displayed X2 — ~X3i X4 = ~~ X5j an d so forth. This 
behavior is quite evident in the numerical solutions for 
S w 1CT 3 . 

Parenthetically, we should note that there is some sim- 
ilarity between the problem that has been solved in this 
section, and a class of problems that has arisen in the con- 
text of magnetic multilayer systems, namely the determi- 
nation of energy (or free energy) minimizing spin confix, 
urations which have variation only in one dimension.Q'El 
Our problem is actually much simpler. By restricting our 
attention to the linear susceptibility the calculation re- 
duces to solution of a linear system, which allows for a 
straightforward, if messy, analytic solution. It is also the 
case that for magnetic multilayers the interest lies in sys- 
tems with large "single-ion" (where an ion corresponds to 
a magnetic layer) XY anisotropy and longitudinal mag- 
netic fields, so there has been no overlap with the results 
of the present work. 

III. QUANTUM SPINS CHAINS 

Let us now return to S = 1/2 quantum spin chains, 
and consider whether the alternation in the T = local 
susceptibility of the semi-infinite Heisenberg chains per- 
sists when the rotational symmetry of the interaction in 
broken. That is, we will consider again the Hamiltonian 
(|J) but reinterpret the S„ as quantum spins rather than 
classical spins. For the XY model (J 2 = 0), it is possi- 
ble to exactly calculate Xn Z ■ An exact calculation of Xn X 
is trivially possible for the Ising model, and furthermore 
one can construct perturbation expansions for the Xn X m 
powers of A. We will find a susceptibility alternation for 
the XY semi-infinite chain, albeit not so strong as for 
the Heisenberg model. For A < 1 we will find behavior 
reminiscent of the classical Heisenberg- Ising semi-infinite 
chain, except that as A — > 1 the local susceptibilities 
diverge classically but we know they are finite quantum- 
mechanically. 
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A. XY Model 



For open chains it is straightforward to apply the 
Jordan- Wigner transformationO to transform the XY 
model with a field h = hi into a model of noninteracting 
spinless fermions, for which the corresponding Hamilto- 
nian has the form 

JV-l N 

n = J±J2 (4^+1 + 4+i*) - - 5) • ( 18 ) 

i=l i=l 

The ground state is constructed by filling all of the neg- 
ative energy single-particle states. The single-particle 
eigenvalues have the form X m — — h + 2J± cosq m where 
q m = Km/ (N + 1) and m runs over the integers 1 . . . N; 
and the normalized eigenfunctions are given by ip m (n) = 
s/2/(N +1) sin q m n. 

The local magnetization is given by 

(S*) = ~+ £ \Mn)\\ (19) 

{m|A m <0} 

but is not uniquely defined if there is an eigenvalue (cor- 
responding to m — mo) which is exactly zero since in 
that case an infinitesimal change in h can cause a finite 
change in (5*), of \ip mo (n)| 2 . 

At h = 0, for finite- length chains the parity of iV 
strongly affects the susceptibilities. For even N, the 
single-particle eigenvalues have a gap of 0(N~ 1 ) around 
zero, and so the zero-field local susceptibilities equal zero. 
The zero-field local magnetizations also vanish. For odd 
N, A(jv+i)/2 = 0; since the corresponding wave vector 
is 7r/2 the susceptibility has (5-function character (with 
strength of 0(N^ 1 )) at h = on odd- numbered sites and 
vanishes on even-numbered sites. (These peculiarities of 
of chains with open boundary conditions were also noted , 
in a slightly different context, by S0rensen and Affleck.Ll) 

However, it turns out that thermodynamic limit is per- 
fectly well behaved despite the disparate characteristics 
of different parity finite-length chains. In the N — * 00 
limit one can replace the sum in ( [l9| ) by an integral, as 
follows: 

1 2 N 

m=(JV+l)(i- 3 ^ r ) 

it/2 

dq sin 2 q . (20) 

tt/2— h/2J± 



(In the sum there should be a greatest-integer function 
acting on the /i-dependent term in the lower limit, how- 
ever, that has no effect when N — * 00.) Differentiation 
with respect to h leads to 

Xl z = (7rJ ± )- 1 sin 2 ^ , (21) 
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that is, Xn Z vanishes on even-numbered sites and takes 
the value (ttJ±)~ 1 on odd-numbered sites. This cer- 
tainly constitutes an alternation in the local susceptibil- 
ity, one which is weaker than in the quantum Heisenberg 
model but, surprisingly, stronger than in the classical XY 
model. 



B. Expansions about the Ising Model 

At present exact results for the T = local trans- 
verse susceptibilities of semi- infinite S — 1/2 chains are 
only available for the Ising and XY models. In order to 
obtain unbiased approximations one might consider car- 
rying out exact diagonalization or quantum Monte Carlo 
calculations on open iV-site chains. Here we will take 
another approach which avoids the need for an explicit 
N — > oo limit, namely, the calculation and extrapolation 
of high-order convergent perturbation expansions about 
the Ising model. Such expansions have proven very infor- 
mative for infinite, homogeneous systems in both one and 
two dimensions;!] for more general discussions of the tech- 
niques used in high-order convergent perturbation expan- 
sions see Ref. 9. We present here the first results obtained 
using such methods for an inhomogeneous system. The 
calculations were carried out by R. R. P. Singh.E3 

Expansions were calculated for J z Xn, with n = 1 . . . 15, 
up to order A 14 (recall A = J±/J z ). Henceforth let us set 
J z = 1 for convenience. The series for \n is identical to 
the series for the bulk \ U P to order A" -2 , and so it is 
more informative to consider the series for A%„ = Xn — X: 
for which the first nonzero term is at order A™ -1 . A 
subset of the calculation series coefficients, sufficient for 
our present purposes, is presented in Table |. 

The most prominent features of the Ax™ series that 
should be evident in the Table are (i) the aforementioned 
vanishing of terms up to order A™~ 2 , and (ii) the sign 
of the leading term alternates between positive, for odd 
n, and negative, for even n. In fact, the same proper- 
ties hold for the semi-infinite classical Heisenberg-Ising 
model, as one can verify by carrying out the Taylor ex- 
pansion in A of the exact classical solution presented in 
Sec. ffj. The leading term in Axn for the classical model is 
(— A/2)™ _1 /2. Note that the coefficients in that case are 
smaller in magnitude than the corresponding coefficients 
for the S = 1/2 system and the difference increases with 
n, so the susceptibility alternation near the Ising limit is 
apparently enhanced by quantum fluctuations. 

Given the reasonably long series available one might 
hope to obtain reliable estimates of the Xn over the whole 
interval [0,1]. However, series extrapolation turns out 
to be rather difficult. Consider first the bulk suscep- 
tibility. Direct Pade approximants appear to be well- 
behaved, and suggest that x f° r the Heisenberg model 
is in the range 0.107 to 0.111. In fact the exact value 
is 0.1013.. O. The Pade approximants are systemati- 
cally in error, albeit modestly, presumably because there 
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is a strong singularity_La \ as A — > (— 1) + and a weak 
singularity as A — > 1 113. It is possible that related sin- 
gularities might exist in all the A\n, although we have 
not considered that problem in detail. A more substan- 
tial barrier to reliable extrapolations based on the A\ n 
series is that, in comparison with the bulk \ series, they 
are quite irregular and the degree of irregularity increases 
with n. Only for A%i do we dare offer an estimate of its 
value at A = 1: a broad interval from 0.7 to 1.3. In 
fact the approximants are clustered at the ends of the 
interval rather than being spread uniformly, so an esti- 
mate presented in the form A ± B would be somewhat 
misleading. 

IV. CONCLUSIONS 

The phenomenon of local transverse-susceptibility al- 
ternation in semi-infinite spin-chains at T — turns out 
to be quite rich. For classical spins, the difference be- 
tween the local susceptibility and the susceptibility of a 
uniform chain decreases with distance from the chain end 
at fixed anisotropy, but diverges as the anisotropy ap- 
proaches zero at any fixed distance. For classical Heisen- 
berg chains at T = the local susceptibilities are all 
divergent and alternating in sign; thus quantum fluctua- 
tions in isotropic 5 — 1/2 chains have the effect of sup- 
pressing susceptibility alternation, most effectively near 
the end of the chain. These findings lead us to sug- 
gest that the susceptibility alternation in semi-infinite 
5=1/2 Heisenberg chainsa is not an essentially quantum 
phenomenon. They might also lead one to conclude that 
quantum fluctuations always suppress the local suscepti- 
bility alternation in semi-infinite chains, but that turns 
out to be incorrect. For both the XY model and near the 
Ising limit, the 5=1/2 chain has stronger alternation 
than the classical chain. 
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FIG. 1. Plot of the local susceptibilities \n for n — 1 
through 6 for the classical Heisenberg-Ising semi-infinite 
chain, as a function of the parameter A = Jx/Jy . The solid 
line is the bulk susceptibility. 
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